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1.  INTRODUCTION 

Most  of  the  work  conducted  during  the  second  research  year 
concentrated  on  extending  and  refining  the  results  found  during 
the  first  year. 

Based  on  up-dating  of  the  literature  survey  on  mechanics  of 
• bonded  joints,  our  data  and  analytical  solutions  were  compared  to 

the  results  of  others,  which  then  raised  certain  questions  neces- 
sitating additional  study.  The  main  objectives  were  to  clarify 
some  dubious  points  and  to  check  the  controversial  information  on 
this  subject. 

Accordingly,  the  work  was  conducted  as  follows: 

— The  finite  element  method  (FEM)  was  used  to  solve  the  prob- 
lem of  stress  distribution  through  the  adhesive  layer  thickness 
and  closer  to  its  edges. 

- The  thermoelastic  behavior  of  symmetrical  and  non-symmetri- 
cal  doubler  configurations  was  studied  with  the  use  of  approximate 
models,  and  analytical  results  were  verified  experimentally. 

- A direct  method  for  measuring  stress  distribution  along  the 
interlaminar  adhesive  layer  (IAL)  was  developed  and  the  analytical 
results  of  the  first  year  were  confirmed  at  a region  very  close  to 
the  IAL  edges. 


- Improved  techniques  for  the  evaluation  of  strength  character- 
istics and  deformational  behavior  of  bonded  adhesives  were  developed 
to  provide  a data  base  for  nonlinear  and  ultimate  stress  analysis. 

- The  effect  of  adhesive  layer  thickness  on  the  above  charac- 
teristics was  studied. 

At  present,  most  of  the  efforts  focus  on  two  major  topics:  an 
analytical  study  of  the  effect  of  the  nonlinear  (inelastic)  IAL 
stress-strain  relation  on  stress  distribution,  and  an  investiga- 
tion of  failure  mechanism  and  strfength  characteristics  of  doublers 
composed  of  aluminum  and  different  composite  adherends.  It  is 
planned  to  continue  these  subjects  during  the  next  research  year. 

The  current  program  and  future  goals  are  illustrated  schemati- 
cally in  Figs.  36  and  37. 


2 . LITERATURE  SURVEY 


An  up-to-date  bibliography  on  the  subjects  of  bonded  joints 
and  interlaminar  behavior  is  provided  in  chronological  order  at 
the  end  of  this  report  in  Appendix  1 which  classifies  this 
literature  according  to  several  categories. 

Most  of  the  literature  reviewed  here  deals  with  joints  composed 
of  rigid  adherends  and  flexible  polymeric  adhesives.  The  adherends 
were  metal-to-metal  (aluminum  and  titanium) , composite-to-composite , 
and,  less  frequently,  metal-to-composite. 

The  mechanical,  analytical  representation  consisted  of  main 
thin  beams  or  plates  on  elastic  substrate  in  bending.  The  common 
structural  models  were  single  and  double  lap  joints  and,  mainly 
in  the  case  of  composite  material  bonding,  scarfed  and  stepped 
joints. 

Linear  elastic  behavior  of  both  adherends  and  adhesives  and 
a state  of  plane  strain  in  the  adhesive  and  of  plane  stress  in 
the  adherends  were  assumed  in  most  cases.  The  common,  approximate 
assumptions  were  uniform  stress  distributions  through  the  adhesive 
thickness,  negligible  transverse  shear  deformations,  and  axial 
normal  interlaminar  stresses.  Consequently,  two  main  variables 
were  assumed  to  have  a predominant  role  in  adhesive  layer  behavior, 
namely,  axial  shear  and  lateral  normal  interlaminar  stresses. 

Boundary  conditions  were  zero  forces  and  moments  at  the  adherend 
ends  and  zero  shear  stress  at  the  adhesive  free  edges  in  the  case 
of  square  geometry.  In  several  closed  form  solutions,  the  last 
condition  was  not  satisfied  because  the  analytical  model  did  not 
have  enough  degrees  of  freedom. 

Earlier  works  neglected  the  presence  of  lateral  normal  stresses 
and  dealt  only  with  axial  shear  along  the  IAL . Volkersen  (1938) , 
in  his  first  study  of  the  problem  [2] , obtained  a simple  expres- 
sion for  the  shear  distribution  of  a single  lap  joint  (SLJ)  com- 
posed of  uniform  adherends  of  different  materials  and  thicknesses. 
According  to  Volkersen,  maximum  shear  stress  is  attained  at  the 
adhesive  free  edge;  the  IAL  stress  level  is  mainly  dependent  on  a 
parameter  consisting  of  the  rigidities  and  thicknesses  of  adherends 
and  adhesive  layers . 

DeBruyen  (1944)  [3]  demonstrated  the  improvement  achieved  in 

stress  distribution  by  varying  the  thickness  of  adherends  down  to 
zero  at  the  edge.  Basing  his  analytical  prediction  on  Volkersen’ s 
solution,  he  found  the  dependence  of  joint  strength  on  the  IAL 
thickness/length  ratio  to  be  the  governing  relationship  for  an 
assessment  of  the  performance  of  bonded  joints. 


Damarkless  (1955)  [8]  extended  Volkersen's  solution  to 

cover  the  effects  of  shear  deformation  in  the  adherend.  An  earlier 
publication  by  Inglis  (1923)  [1]  indicated  the  existence  of  pro- 

nounced transverse  stresses  and  stress  concentration  at  the  joint 
edges.  This  problem  was  studied  much  later  by  Adams  and  Peppiatt 
(1973)  [40]  , who  revealed  the  presence  of  high  transverse  shear 

stresses  o.i  the  IAL  which  were  due  to  the  Poisson  effect  of  the 
adherend. 

The  breakthrough  in  this  field  was  achieved  by  the  classical 
contribution  of  Goland  and  Reissner  (1944)  [4]  , who  provided  a 

closed  form  analytical  solution  for  single  lap  joints  with  ex- 
pressions of  shear  and  lateral  normal  stress  distributions  through 
the  IAL.  They  presented  a model  of  thin  adherend  plates  of  equal 
thickness  and  rigidity  under  cylindrical  bending.  Here  again, 
though,  the  boundary  conditions  of  zero  shear  stress  at  the  IAL 
edges  was  unsatisf ied.  An  approximate  solution  based  on  their 
work  was  given  by  Plantema  (1949)  [5]  ; it  included  the  effect  of 

external  moment  on  IAL  shear  stresses,  but  neglected  this  effect 
on  the  respective  normal  stresses. 

Sherrer  (1957)  [10]  extended  the  work  of  [4]  to  cover  ad- 

herends  of  different  materials  by  assuming  a plane  stress  for  the 
adherends  and  a plane  strain  for  the  adhesive. 

Cornell  (1953)  [7]  presented  a model  of  beams  representing 

the  adherends  not  affected  by  the  adhesive  layer,  which  was 
represented  by  a series  of  tensile  and  torsional  springs. 

Iyengar  and  Alwar  (1961)  [13]  attempted  to  solve  the  model 
of  [4]  by  fulfilling  boundary  conditions  for  the  free  IAL  edges. 
Kunzi  and  Stevens  (1963)  [15]  corrected  the  solution  of  [4]  and 

extended  it  to  cover  adherends  of  different  characteristics. 

Volkersen  (1965)  118],  in  a later  work,  suggested  a solution 

which  included  both  shear  and  normal  stresses  and  satisfied  the 
IAL  free  edge  conditions. 

The  above-mentioned  solutions  following  [4]  were  based  on 
many  approximate  assumptions  and  over-simplified  models.  This 
was  due  to  the  complexity  of  an  exact  model.  It  was  later  found 
that  a few  assumptions  led  to  an  incorrect  solution  at  the  cri- 
tical region  near  the  edges.  Numerical  solutions  were  suggested 
in  order  to  obtain  a more  exact  solution,  and  the  major  effort 
was  devoted  to  applying  the  finite  element  method  (FEM)  to  solve 
the  problem.  At  first,  this  method  was  applied  for  checking  ex- 
perimental tests  to  determine  joint  strength  and,  especially,  the 
effect  of  edge  geometry.  Later,  the  whole  IAL  stress  distribution 
was  derived  by  the  FEM. 


In  recent  years,  several  works  were  published  which  used 
FEM  to  solve  the  problems  of  SLJ  and  DLJ  models.  Among  the 
first  was  the  report  by  Kutscha  et  al.  (1969)  [25],  which  is 

based  on  available  FEM  programs.  The  solution  derived  is  for 
states  of  plane  stress  and  plane  strain.  The  model  used  was  the 
same  as  those  of  [2]  and  [4] , and  the  results  provided  a measure 
for  the  accuracy  of  the  closed  form  analytical  solution. 

Wooley  and  Carver  (1971)  [33]  developed  a general  solution 

for  SLJ,  based  on  FEM,  which  agreed  well  with  [4]  but  did  not 
account  for  the  free  edge  boundary  conditions . 

Dickson  et  al . (1972)  [36]  derived  analytical  solution  for 

joints  with  orthotropic  adherends;  they  also  developed  a computer 
program  based  on  FEM  which  enables  one  to  deal  with  different 
joint  geometries  and  non-linear  IAL  stress-strain  behavior. 

A similar  attempt  to  deal  with  non-linear  behavior  was 
carried  out  by  Grimes  et  al.  (1972)  [39],  who  suggested  a 

numerical  solution  to  the  problem  of  joints  composed  of  orthotropic 
FRP  adherends . 

Contributions  by  Barker  and  Hatt  (19  73)  [47]  , and  Adams  and 

Peppiatt  (1974)  [51] , dealt  with  the  problem  of  singularity  at 

the  IAL  edges  by  means  of  the  finite  element  method.  At  the 
same  time,  works  by  Hart-Smith  (1973)  [42-46]  and  Renton  and 

Vinson  (1974)  (41,58)  suggested  analytical  solutions  to  the 

problem  of  structural  bonded  joint  including  anisotropic  compo- 
site adherends. 

An  attempt  to  tackle  the  problem  of  stress  distribution 
through  the  IAL  thickness  was  undertaken  by  Pirvics  (1974)  [5] 

using  finite  difference  methods.  His  solution  may  violate  the 
common  assumption  of  uniform  shear  and  tension  through  the  IAL 
that  was  postulated  by  many  authors. 

In  recent  years,  much  more  attention  has  been  given  to  the 
present  subject;  the  scope  was  expanded,  and  efforts  were  more 
orinted  towards  the  applied  aspects  of  analysis  and  the  design 
of  structural  bonding. 

To  sum  up,  the  numerous  contributions  on  the  subject  during 
the  past  decade  may  be  classified  into  the  following  topics: 

I.  Mechanics  of  joints  composed  of  orthotropic  composite  ad- 
herends [reference  numbers:  20,  23,  25,  29,  36,  37,  39,  41, 

42,  43,  55,  56,  58,  59,  63] . 

II.  Fracture  mechanics  approach  and  testing  methods  [reference 

numbers:  49,  53,  60,  64]. 
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III.  Environmental  effects  on  bonded  joint  performance  (includ- 
ing thermoelastic  behavior),  [reference  numbers:  21,  35, 

49,  50,  54,  64]. 

IV.  Dynamic  behavior  (fatigue)  of  structural  joints  [reference 
numbers:  21,  32,  36,  41,  54,  55,  56,  63]. 

V.  Effect  of  non-elastic  adhesive  characteristics  on  bonded 
joint  performance  [reference  numbers:  29,  37,  38,  39,  42, 

46,  631. 

VI.  Effect  of  varying  geometry  and  adhesive  composition  at 
the  IAL  boundaries  on  structural  performance  of  the  joint, 
(this  may  provide  guidelines  for  optimization  of  structural 
bonding)  [reference  numbers:  17,  3 , 40,  47,  51,  611. 

VII.  Analysis  of  joints  (stepped  and  scarfed)  with  complex  con- 
figurations to  fit  design  requirements  [reference  numbers: 

34 , 45,  46,  47,  55,  59 , 61] . 

VIII.  Design  considerations  and  evaluation  of  design  allowables 
for  structural  bonded  joints  [reference  numbers:  19,  23, 

24,  29,  31,  37,  46,  54,  58]  . 

The  last  topic,  VIII,  is  based  on  the  work  by  Hughes  and 
Rutherford  (1968)  [24]  , who  developed  a testing  methodology  for 

direct  measurement  of  stress-strain  behavior,  modulus  and  strength 
of  the  adhesive  layer  in  situ.  In  the  course  of  their  investiga- 
tion, they  evaluated  the  mechanical  characteristics  in  shear  and 
tension  of  different  structural  adhesives  of  varying  thicknesses. 

Numerous  articles  are  available  on  the  chemical  and  physical 
aspects  of  the  problem,  as  well  as  on  the  more  technological  as- 
pects involved  in  the  fabrication  and  inspection  of  adhesive 
joints.  These  contributions  are  beyond  the  scope  of  the  present 
survey,  which  is  limited  to  the  subject  of  the  mechanical  behavior 
of  structural  bonded  systems. 


3.  TWO  DIMENSIONAL  INTERLAMINAR  STRESS  DISTRIBUTION  IN  THE 
ADHESIVE  LAYER  OF  A SYMMETRICAL  DOUBLER  MODEL. 


3.1  Introduction 

In  order  to  derive  a closed  form  analytical  solution  for  stress 
distribution  of  the  IAL , the  previous  report  [65]  postulated  cer- 
tain simplified  assumptions,  which  were  shared  by  other  authors  [4, 
7,15,18,41]:  uniform  shear  and  normal  stress  distribution  through 
the  thickness  direction  (z)  . 

The  neglect  of  axial  normal  stress  (o^)  in  the  adhesive  is  also 
common  to  many  publications  on  bonded  lap  joints.  These  simpli- 
fied assumptions  and  the  fact  that  the  analytical  model  does  not 
have  sufficient  degree  of  freedom  lead  to  several  inconsistencies 
in  equilibrium  relationships  and  violate  the  boundary  condition 
of  zero  shear  stress  at  the  free  IAL  square  edges. 

Because  of  the  mathematical  complexity  required  of  a more 
rigorous  solution,  and  the  fact  that  a realistic  bond-line  is 
very  thin  relative  to  adherend  thickness,  only  a few  studies  have 
been  made  on  IAL  lateral  stress  distribution.  Pirvics  [52  ob- 
tained a solution  by  finite  difference  minimization  of  the  internal 
energy  distribution  for  typical  simple  lap  and  butt  joint  models. 
His  solution,  though  suffering  from  inconsistencies,  clearly  in- 
dicates the  pronounced  non-uniformity  in  the  distribution  of  both 
shear  and  normal  stresses  through  the  adhesive  thickness. 

The  present  study  is  intended  to  deal  with  the  problem  by 
means  of  the  finite  element  method.  The  main  purposes  were  to 
check  analytical  solutions  and  to  provide  basic  data  for  an  in- 
vestigation of  strength  and  failure,  for  which  the  local  princi- 
pal stress  components,  rather  than  their  average  values,  are 
required. 


3.2  Model  Representation 

The  model  for  the  finite  element  solution  of  the  SMD  is 
based  on  orthogonal  and  triangular  elements  with  a uniform  strain 
field.  The  elements  satisfy  the  compatability  conditions  at 
each  of  the  material  points,  as  well  as  equilibrium  conditions  in 
each  element.  The  mechanical  behavior  is  limited  to  the  elastic- 
linear  range  and  to  the  two-dimens icr.al  case.  A state  of  plane 
strain  was  assumed  for  the  IAL  (x-z  plane)  because  of  its  high 
width  to  thickness  ratio  and  its  low  adhesive-to-adherend  stiff- 
ness ratio.  Although  this  assumption  is  less  accurate  when  re- 
lating to  the  adherend  layers,  it  is  justified  for  the  present 
study,  which  focuses  on  the  IAL.  The  results  of  preliminary  com- 
parative studies  between  plane  stress  and  plane  strain  conditions 
for  the  adherends  exhibited  only  small  differences. 
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IAL  behavior  is  best  represented  by  the  SMD  structural 
model  shown  in  Fig.  1.  This  model  was  also  selected  for  the 
FEM  program  because  of  the  following  advantages: 

- Double  symmetry  in  geometry,  materials,  boundary,  and 
loading  conditions,  which  allows  treatment  of  one  quarter  of  the 
model  only  (in  contrast  with  the  SLJ  model) . 

- Wide  separation  between  external  loading  region  and  the 
critical  IAL  zones  where  maximum  stresses  and  strains  are  located. 
This  allows  the  isolation  of  the  IAL  boundary  zone  which  is  of 
main  interest  for  both  experimental  and  analytical  study. 


- Continuity  of  stress  and  strain  fields  over  all  the  IAL 
region,  which  reduces  the  number  of  boundary  conditions  (in  con- 
trast to  the  case  of  the  SLJ  model) . 


The  model  for  the  present  specific  case  is  composed  of 
aluminum  adherends  and  epoxy  adhesive,  the  characteristics  of 
which  are  given  in  Table  1 (test  model  A21) . 

The  FEM  program  consisted  of  three  steps.  In  the  first 
step,  uniform  stress  distributions  through  IAL  thickness  was  as- 
sumed and  the  IAL  was  divided  into  rectangular  elements  having 
the  IAL  lateral  dimensions  (Fig.  2a> . 

In  the  second  step,  the  IAL  thickness  was  divided  into  two 
longitudinal  strips  (Fig,  2b) . in  the  final  step,  which  aimed  at 
deriving  the  stress  distribution  through  the  IAL  thickness,  it 
was  divided  into  4 longitudinal  strips  as  shown  in  Fig.  2c. 


3.3  FEM  vs.  Analytical  Solution 


The  FEM  was  first  examined  by  comparing  its  results  (based 
on  the  model  of  Fig.  2A)  to  the  analytical  closed-form  solution 
of  ref.  [65].  Good  agreement  was  found  between  numerical  and 
analytical  results  both  for  the  normal  stress  distribution  (Fig. 
3)  and  for  the  shear  stress  distribution,  except  at  the  region 
close  to  the  IAL  edges  (Fig.  4).*)  A preliminary  study  of  the  ef- 
fect of  a finer  network  m the  axial  direction  showed  only  slight 
differences  in  results. 


The  second  step  revealed  different  stress  patterns  close  to 
the  edges.  Shear  stress  distribution  was  found  to  reverse  its 
slope  at  the  boundary  zone  (X^O.98)  and  to  drop  to  zero  towards 
the  IAL  edges  (Fig,  3).  Normal  stress,  on  the  other  hand,  in- 
creased steeply  towards  the  edge  of  the  IAL  layer  close  to  the 


*)  In  all  stress  distributions  (Figs.  3-17)  the  stresses  are  nor- 
malized with  respect  to  average  applied  stress  acting  at  the  cen- 
tral adherend  pc=F/2h  . Accordingly:  T XZi=i X2i/Pc • BXi=0xi/pc 
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central  adherend  but  reversed  their  slope  and  dropped  to  zero 
towards  the  free  edges  of  the  IAL  layers  close  the  external  adher- 
ends  (Fig.  4) . 

3.4  Two-Dimensional  Stress  Distribution 

The  third  step  of  the  FEM  solution,  in  which  the  IAL  was 
divided  into  4 sub-layers,  yielded  a two-dimensional  stress  dis- 
tribution for  both  adhesive  and  adherend.*) 

3.4.1  Stress  Distribution  in  Adherends 

The  predominant  stress  component  obtained  at  the  adherends 
is  the  axial  normal  one  - (Figs.  5-6)  , which  was  higher  by 

an  order  of  magnitude  than  the  other  adherend  stresses,  a . 

and  t . (Figs.  7-10);  a . was  almost  uniform  through  the  adfierend 

■ i * i XZ 1 XI 

thickness . 

In  the  case  of  metallic  adherends,  in  which  other  stresses 
besides  a . are  negligible,  the  assumption  of  uniform  stress  dis- 
tribution through  the  adherend  thickness  was  justified.  However, 
with  adherends  composed  of  multilayer,  composite  laminates,  in 
which  interlaminar  stiffness  and  strength  characteristics  are 
significantly  lower  than  axial  ones,  adherend  lateral  stresses 
t and  a . have  a major  role  and  cannot  be  neglected.  These 

XZ  X Z 1 

stresses  exhibited  pronounced  variations  through  the  thickness. 

In  this  case,  a finer  mesh  for  the  FEM  solution  is  essential,  es- 
pecially if  strength  characteristics  are  required. 

The  axial  distribution  of  axial  normal  stresses  through  the 
external  layer  of  adherend  2 exhibited  the  pattern  which  was 
found  both  analytically  and  experimentally  in  ref.  [65]  (Fig.  30), 
namely  a change  from  tension  to  compression  close  to  the  adherend 
edge  (Fig.  5) . 

3.4.2  Stress  Distributions  Through  Adhesive 

In  the  IAL,  all  stress  components  are  of  the  same  order  of 
magnitude  and  exhibited  a pronounced  variation  through  the  thick- 
ness (along  the  z axis);  (Figs.  11-13,  14-16).  The  axial  distri- 
bution of  shear  stress  (t  ) showed  a trend  similar  to  that  in 
Step  2,  i.e.,  slope  reversal  at  the  'boundary  zone'  (from  each 
edge  up  to  a distance  of  1%  of  the  IAL  length) , and  a drop  towards 
zero  at  the  free  IAL  edges.  This  finding  is  consistent  with 
boundary  conditions  but  contradicts  the  closed-form  analytical 
solution  (Figs.  12,15). 

Axial  distribution  of  lateral  normal  stresses  (a  ) along 
the 'boundary  zone'  showed  a pronounced  divergence  from  their 
_ 

The  results  obtained  in  the  present  report  were  related  to  a 
model  of  isotropic  adherends.  The  FEM  program,  which  is  given  in 
Appendix  2,  is,  however,  adequate  for  orthotropic  adherend  cases 
as  well. 


9 


average  reference  distribution,  whereas  the  normal  stress  through 
the  external  IAL  layer  rose  steeply  towards  the  edge  far  beyond 
its  slope  and  dropped  towards  zero  at  the  free  edge  (Figs.  13,16). 

Axial  distribution  of  axial  normal  stresses  (axo)  was  uniform 
throughout  the  "middle  zone"  (98%  of  the  IAL  length)  but  dropped 
steeply  at  the  "boundary  zone"  towards  zero  at  the  free  edges 
(Fig.s  11,  14) . 

It  may  be  concluded  that  axial  interlaminar  stress  distribu- 
tions (ozo,  txzo)  derived  by  FEM  are  in  close  agreement  with  the 
closed-form  solutions  along  the  "middle  zone".  The  major  diver- 
gence between  the  ftumerical  and  closed-form  solutions  resulted 
from  the  assumption  of  uniform  lateral  stress  distributions. 

As  shown  in  Figs.  15,16,  such  an  assumption  did  not  hold  up  at 
the  "boundary  zone"  in  cases  of  shear  and  normal  lateral  stresses. 
The  neglect  of  axial  normal  stresses  through  the  IAL  is  also  un- 
justified as  these  stresses  attained  a level  having  the  same 
order  of  magnitude  as  other  stress  components  (Fig.  14) . 


Lateral  normal  stress  distribution  at  the  boundary  zone  shows 
the  maximum  deviation  from  uniformity  (Figs.  13,16),  being  almost 
zero  close  to  the  external  adherend  and  increasing  steeply  an4 
linear! ly  to  its  maximum  level  towards  the  central  adherend. 
Similarly,  both  IAL  shear  and  axial  normal  stresses  attained  their 
maximum  value  close  to  the  intersection  of  the  free  IAL  edge  ^nd 
the  upper  surface  of  the  central  adherend  (point  m of  Fig.  1) . 

3.5  Discussion 

Given  the  values  of  the  three  stress  components  at  the  IAL 
boundary  zone,  their  principal  value  can  be  derived.  Fig.  17 
shows  the  axial  and  lateral  distributions  of  the  principal  stresses 
emphasizing  the  major  trends  found  for  the  stress  components. 

The  expected  finding  of  the  location  of  maximum  stress  being  close 
to  point  m is  also  demonstrated. 

In  the  case  of  brittle  IAL  characteristics,  the  data  given  by 
the  principal  stress  distribution  may  provide  the  basic  clue  for 
failure  mechanism,  the  delamination  process,  and  strength  evalua- 
tion for  the  overall  bonded  structure  as  represented  by  the  SMD 
model.  The  finite  element  method  seems  to  provide  an  accurate 
solution  which  conforms  well  to  basic  mechanical  and  boundary 
conditions.  The  interrelation  among  the  different  functions 
found  for  the  different  stresses  can  be  examined  by  satisfying  the 
classical  equilibrium  conditions. 

According  to  the  first  equilibrium  equation 


(1) 


.«* 
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which  means  that  txz  is  uniformly  distributed  along  the  z axis 
in  the  region  where  ax  is  uniformly  distributed  along  x.  This 
interrelation  prevails  at  the  "middle  zone"  of  the  IAL,  as  is 
shown  when  Fig.  15  is  compared  with  Fig.  14. 

According  to  the  second  equilibrium  equation. 


do 


z 


cTz 


dx 


xz 

dx 


(2) 


which  means  that  where  the  TX2  distribution  along  x reverses  its 
slope,  the  a function  along  z will  also  reverse  but  in  the  oppo- 
site direction.  This  change  occurs  simultaneously  for  the  two 
functions  at  about  X=0.99  within  the  "boundary  zone",  as  demon- 
strated by  a comparison  of  Fig.  15  with  Fig.  16. 

At  the  middle  zone,  the  txz(x)  functions  coincided  (Fig.  15), 
thus  oz  functions  are  supposed  to  be  linear  functions  of  z,  as 
is  clearly  shown  in  Fig,,  16. 

The  above  findings  referred  to  an  idealized  case  of  thick 
adhesive  relative  to  its  adherends.  It  is  expected  that  similar 
distribution  patterns  will  be  found,  though  to  a lesser  degree, 
for  the  more  practical  cases  of  thinner  and  stiffer  adhesives. 

The  most  important  finding  of  the  locations  of  the  maximum 
principal  stresses  is  also  expected  to  be  generally  valid. 

A more  detailed  solution  tor  interlaminar  adherend  stresses 
is  essential  for  the  strength  analysis  of  orthotropic  composite 
bonded  systems  - 

In  order  to  obtain  more  realistic  solutions  which  provide 
quantitative  data  for  design  allowables,  the  effects  of  nonlinearity 
in  an  adhesive  stress-strain  relationship  and  the  effect  of  edge 
geometry  have  to  be  studied. 

Some  of  the  above  topics  are  currently  being  treated  by  the 
finite  element  method  and  wiil  be  presented  in  the  future. 
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3.6  Conclusions 


The  following  conclusions  emanate  from  the  solutions  ob- 
tained by  applying  the  finite  element  method  to  the  symmetrical 
doubler  model,  composed  of  aluminum  adherends  bonded  by  thick 
adhesive  layers  (IAL) . Most  of  the  conclusions  concern  the 
stress  distributions  through  the  IAL,  which  is  the  main  subject 
of  the  present  investigation. 

1.  The  distribution  of  shear  and  lateral  normal  stress  through 
the  IAL  thickness  was  approximately  uniform  within  the  "middle 
zone"  (-0 . 98<x<0 . 98) . 

2.  Axial  normal  stress  distribution  [a  (x) ] was  uniform  within 
the "middle  zone"  but  varied  linearly  through  the  IAL  thickness. 

3.  For  all  interlaminar  stress  functions,  drastic  variations 
occurred  at  the  "boundary  zone". 

4.  The  distribution  of  IAL  axial  shear  stress  attained  its 
maxima  within  the  "boundary  zone",  reversed  its  slope  and  dropped 
towards  zero  at  the  free  edge.  A similar  trend  was  found  for  oz 
distributions  close  to  the  external  adherends . 

5.  The  normal  lateral  stress  function  a (z)  exhibited  an  approxi- 
mately linear  distribution,  which  reversed  its  slope  abruptly 
within  the  "boundary  zone". 

6.  Maximum  values  of  all  stresses  and  their  principal  counter- 
parts were  attained  at  the  "boundary  zone",  close  to  the  central 
adherend  interface. 

7.  The  different  patterns  which  characterized  stress  functions 
along  the  z and  x axes  were  consistent  with  basic  differential 
equilibrium  equations  and  boundary  conditions  of  the  IAL. 

8.  The  average  interlaminar  stress  levels  were  significantly 
lower  than  their  actual  individual  extreme  levels;  thus,  they 
could  not  be  used  to  predict  the  doubler  strength. 

9.  Because  of  the  narrow  dimensions  of  the  boundary  zone,  its 
stress  behavior  will  not  be  reflected  in  the  overly  deformational 
behavior  of  the  bonded  structure,  nor  could  this  behavior  be 
directly  detected  experimentally.  It  is,  however,  predominantly 
influential  on  the  ultimate  behavior  of  the  IAL  and  has  to  be 
considered  when  failure  mechanism  and  strength  characteristics 
are  discussed. 


12 


wm 


t- 

V 

! ' *■ 
; i 

! 


I . '4 


4.  DIRECT  DETERMINATION  OF  AVERAGE  INTERLAMINAR  STRESS 
DISTRIBUTIONS  IN  A POLYMERIC  ADHESIVE  LAYER. 


4 . 1 Introduction 


The  real  stress  distribution  through  IAL  is  still  undecided 
because  of  the  controversy  over  the  different  solutions  of  the 
problem.  Unfortunately,  there  is  no  satisfactory  experimental 
means  of  direct  characterization  of  these  stresses.  Matting  116] 
and  Kreiger  [49]  measured  directly  only  the  shear  displacement  of 
SLJ  and  compared  it  with  solutions  of  [2]  and  [4] . Most  attempts 
at  such  characterization  employed  indirect  measurement  of  strains 
at  the  external  adherend  surfaces  by  means  of  strain  gages  [65,42] 
or  photoelastic  [65,36,9],  or  Moire  [70]  methods.  In  all  cases, 
the  patterns  were  obscured  by  the  interference  of  other  factors 
such  as  stress  concentrations,  edge  effects  of  the  photo-stress 
coating,  etc.  Moreover,  the  interlaminar  stress  distribution 
cannot  be  conclusively  established  from  data  on  adherend  strains. 
This  problem  is  crucial,  considering  that  the  location  of  the 
joint,  where  extreme  stresses  develop  and  enhance  delamination, 
is  very  close  to  the  edges  of  the  adhesive  layer.  The  object  of 
the  present  investigation  is,  therefore,  twofold: 

a.  To  develop  a method  for  direct  determination  of  the 
interlaminar  strain  and  stress  distribution  which  will 
provide  a check  on  the  validity. of  existing  analytical 
solutions ; 

b.  to  devise  a tool  for  the  direct  measurement  of  the 
critical  stresses  which  control  the  ultimate  performance 
of  the  joint  or  laminate  under  different  loads  and  envi- 
ronmental conditions. 

4.2  Analytical  Background 

A simplified  expression  of  IAL  shear  stress  distribution 
for  the  SMD  model  (Fig.  1)  is  given  in  eq.  (3) : 


T*  - h<  <P  sinh<fr~x 
c coshc(> 


(3) 


where  r* 


T 

Pu 


is  the  normalized  shear  stress; 


<t>* 


C*Gp 
h o h 3 E i 


(1  + M) 
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y 
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and 


Po  = ikitt+ut  (6) 

is  the  reference  axial  normal  stress  (acting  at  the  external 
adherend  layer  -(EAL)  midpoint  under  the  action  of  F) . The 
derivation  of  eq.  (3)  and  the  more  exact  solution  for  shear 
stress  distribution  is  given  in  ref.  [65], 

The  similarity  of  the  simplified  to  the  more  exact  solution  for 
a case  of  multimaterial  doubler  model  is  shown  in  Fig.  9 of 
ref.  [65]. 

The  derivation  of  the  expression  for  interlaminar  lateral 
normal  stresses,  a (°),  is  more  complex  than  that  for  shear  stress 
analysis.  It  involves  the  solution  of  differential  equations  of 
the  6th  order,  as  is  fully  described  in  ref.  [65].  Typical  dis- 
tribution curves  of  a * for  different  boundary  conditions  are 
shown  in  Fig.  10  of  this  reference.  (a**az  (°  Vpo  is  the  normalized 
expression  for  IAL  lateral  stresses.) 

4 . 3 Experimental  Procedure 

Shear  «train  y and  shear  stress  t are  obtainable  from  the 
shear  displacement  Au,  i.e. , the  differential  displacement  of  the 
adherend  surfaces  (Fig.  19): 


Y 


(7) 


Similarly,  by  measuring  the  net  lateral  displacement,  normal 
lateral  stresses  are  obtainable  from  eq.  (8); 


CTo  = {yj-  E0  (8) 

Instron  electromechanical  extensometers  were  modified  and  used 
for  measuring  Au  and  wo , respectively. 

Measuring  Devices . The  shear-displacement  extensometer 
shown  in  Fig.  20  is  a modification  of  the  Instron  Extensometer 
Type  G51-11 . It  consisted  of  two  parts,  each  attached  to  a point 
at  the  mid-thickness  of  the  adherends,  all  of  which  were  origin- 
ally located  on  the  same  line  (along  the  z-axis) , but  which  were 
shifted  relative  to  each  other  under  axial  loading,  thereby 
creating  the  shear  displacement.  The  normal  displacement  distri- 
bution was  obtained  in  the  same  manner,  utilizing  the  normal- 
displacement  extensometer  shown  in  Fig.  21,  which  is  a modifica- 
tion of  Instron  Extensometer  G57-11  (originally  used  for  measuring 
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the  Poisson  ratio) . 

The  displacements  were  measured  within  an  accuracy  of  p.2ym 
The  distributions  were  obtained  by  applying  the  extensometers  at 
different  points  along  the  model  axis. 

Preparation  of  Test  Model.  The  measuring  systems  were 
checked  on  test  models  consisting  of  two  aluminum  adherends  and 
an  epoxy  IAL  having  the  characteristics  given  in  Table  1. 

The  surfaces  of  the  adherends  were  sand  blasted,  cleaned 
with  MEK  and  coated  with  a thin  layer  of  resin,  as  described  in 
detail  in  ref.  [65].  Immediately  following  the  coating,  the 
adherend  strips  were  inserted  into  a special  fixture  which 
ensured  the  exact  spacing  for  the  IAL  phase.  The  external  ad- 
herends were  longer  than  the  adhesive  bond-line  in  order  to 
enable  measurement  of  displacements  at  IAL  edges  (see  Fig.  1) . 

The  epoxy- vers amid  mix  (2:1)  was  prepared  under  vacuum  and 
poured  into  the  gap;  its  temperature  was  kept  at  40 °C  to  reduce 
viscosity  and  ensure  good  wetting.  The  model  was  unmolded  after 
24  hours  and  exposed  for  8 hours  at  40 °C  to  post-cure  the  resin. 

Loading  Procedure . The  model  was,loaded  uniaxially  on  an 
Instron  Tester  at  the  rate  of  about  10  J (min  i) , up  to  failure. 
Shear  displacement  values  were  recorded  continuously.  The  speci 
men  was  unloaded  at  different  loading  levels  as  a check  on  the 
reversibility  of  the  response. 

4.4  Test  Results  and  Discussion 

The  relationship  between  the  shear  displacement  and  the 
reference  axial  stress,  as  shown  in  Fig.  22,  is  characterized 
by  linear  behavior  up  to  approximately  80%  of  the  loading  range. 
Beyond  this  level,  nonlinearity  prevailed  to  the  delamination 
stage . 

The  shear  stresses  at  the  various  points  along. the  IAL  were 
calculated  from  the  measured  shear-displacement  data  with  the 
aid  of  eq.  (7).  The  experimental  interlaminar  shear  stress  dis- 
tribution, shown  in  Fig.  22,  was  compared  with  the  approximate 
analytical  solution  of  eq.  (3).  Considering  the  approximations 
involved,  agreement  between  the  experimental  data  and  the  ana- 
lytical prediction  is  good,  especially  with  regard  to  the  shape 
of  the  distribution  curve  close  to  the  boundary. 

The  tendency  of  the  shear  stresses  to  increase  to  their 
maximum  value  at  the  IAL  edge  was  in  agreement  with  the  analyti- 
cal solutions  of  ref.  [4]  for  bonded  joints  and  of  ref.  [65]  for 
doubler  models.  It  conflicted,  however,  with  the  boundary  condi 
tion  of  zero  shear  stress  at  the  free  IAL  edqe  and  with  solu- 
tions based  on  the  finite-element  method  (Fig.  15) . 
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In  order  to  obtain  the  exact  value  of  lateral  displacement 
resulting  from  interlaminar  normal  stress  o (°),  the  contribution 
of  Poisson's  effect  to  the  overall  displacement  must  be  determined. 
The  Poisson  displacement  data  hardly  vary  along  the  model  axis  and 
may  be  calculated  or  measured  separately. 

The  relationship  between  lateral  normal  displacement  and 
the  reference  axial  stress  is  shown  in  Fig.  22,  where  similar 
trends  to  those  of  the  shear-displacement  data  are  apparent,  i.e., 
linearity  within  80%  of  the  loading  range. 

The  lateral  normal  stresses  at  the  various  points  along  the 
IAL  were  calculated  from  the  net  lateral  displacement  data  with 
the  aid  of  eq.  (8) . 

The  experimental  interlaminar  normal  stress  distribution 
(Fig.  24)  was  compared  to  the  analytical  solution  given  in  eq.  (33) 
of  ref.  [65]  . The  agreement  between  the  experimental  data  and  the 
analytical  prediction  was  good,  especially  with  regard  to  the 
shape  of  the  distribution  curves  close  to  the  IAL  edges. 

4.5  Conclusions 


A method  of  direct  measurement  of  shear  and  normal  displace- 
ments and  a simple  way  of  determining  strain  and  stress  distribu- 
tions along  an  interlaminar  adhesive  layer  were  examined  success- 
fully on  a symmetrical  doubler  model. 

The  experimental  interlaminar  stress  data  were  in  agreement 
with  available  analytical  solutions  for  this  model.  The  advan- 
tages of  the  direct  method  compared  with  others  are  as  follows: 

a.  It  is  based  on  direct,  local  measurement  of  the  deforma- 
tions of  the  adhesive  layer  (as  against  indirect  measure- 
ments of  adherend  strains) . 

b.  It  is  a simple  method  utilizing  the  standard  accessories 
of  an  Instron  extensometer  with  only  slight  modifica- 
tions . 

c.  The  accuracy  in  measuring  displacements  to  within  0.2y  2 

permits  detection  of  stress  variations  of  about  0.2  kg/mni  , 
which  is  less  than  5%  of  the  ultimate  shear  and  stress 
delamination  level. 

d.  The  continuous  recording  of  interlaminar  adhesive  strains 
against  axial  load  permits  monitoring  of  IAL  behavior 
from  the  linear  elastic  range  to  the  nonlinear  inelastic 
range  up  to  initiation  of  the  failure  process. 

This  study  may  reveal  the  actual  mechanism  of  adhesive  and 
cohesive  failure  modes  which  are  eventually  responsible  for  the 
delamination  of  laminates,  joints,  and  other  bonded  systems  that 
are  subjected  to  mechanical  and  thermal  loading. 
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5.  THERMOELASTIC  BEHAVIOR  OF  MULTIMATERIAL  DOUBLER  MODELS 

The  analytical  expressions  in  Ref  [65],  (eqs.  63,  68,  43), 
covered  the  effect  of  thermal  changes  on  mechanical  behavior 
with  respect  to  elastic  and  thermoelastic  parameters. 

The  solution  was  based  on  the  assumption  that  the  coefficient 
of  thermal  expansion,  and  of  other  material  parameters,  were  not 
affected  by  thermal  changes.  The  solution  also  neglected  trans- 
verse behavior  and  its  interrelation  with  axial  behavior.  Within 
these  limitations,  which  apply  for  an  isothermal  regime,  there 
is  a wide  range  of  materials,  geometries,  and  thermoelastic 
parameters  for  which  the  solution  is  valid. 

The  thermoelastic  behavior  of  symmetrical  and  nonsymmetrical 
doubler  models  (SMD  and  NMD,  respectively)  were  studied  with  a 
view  to  checking  the  validity  of  analytical  formulations  of  the 
stress  distribution  by  means  of  strain-  and  deflection-measurements 
along  the  model's  external  facings. 

5.1  Thermoelastic  Behavior  of  a Symmetrical  Doubler  Model  (SMD) 


5.1.1  Analytical  Background. 


The  normalized  approximate  expression  for  the  external  strain 
distribution  of  the  EAL  is  given  by  eq.  9 (eq.  (73)  of  ref.  65) : 

e*  = 1 - C-°-Si  e2  (9) 

cosh<J)  Pc 


where 


E,h, (a, -as ) T 
ha (1+u) 


For  non-loaded  specimens,  eqs.  (9)  and  (10)  yield  the  strain 
distribution  along  the  EAL  external  facing: 

11  - Hi?1  (U> 

E h 

where  u = — v ; a,,  a*,  EJf  E2,  are  the  thermal  coefficients  and 
r,ina 

the  elastic  moduli  in  the  x-direction,  respectively. 

Eqs  9-11  are  based  on  the  assumption  of  negligible  transverse 
stresses,  i.e.: 

o < < o 

zy  zx 


t 


In  the  thermoelastic  case,  such  conditions  require  near-uniformity 
of  the  thermal  expansion  coefficient  in  the  transverse  direction, 
i.e.,  ctly  « u2y. 

5.1.2  Experimental  Procedure. 


Model  A02,  which  had  served  previously  for  experimental 
verification  of  eq.  (9)  in  the  the  case  of  mechanical  axial  load- 
ing, was  also  used  in  the  thermal  tests.  Its  basic  characteris- 
tics are  given  in  Table  2. 


The  axial  thermal  coefficients  ax  and  012  were  measured  by 
strain  gage,  whereas  a is  given  by  xthe  manufacturer  for  the 
GRP  1002  prepregs  usedY1  It  can  be  seen  from  Table  2 that  the 
near-uniformity  requirement  was  satisfied.  The  test  was  carried 
out  as  follows: 


Fourteen  gages  were  glued  to  the  aluminum  external  EAL  sur- 
faces of  the  test  model,  as  shown  in  Fig.  25.  The  specimen  was 
exposed  to  changes  in  temperature  within  a thermal  chamber  that 
ranged  from  30 °C  to  60 °C.  (Higher  temperatures  were  expected 
to  affect  the  elastic  parameters  of  the  polymers.)  Strains  were 
recorded  simultaneously  at  the  different  gage  points.  The  tem- 
perature was  raised  and  lowered  several  times  to  check  reproduci- 
bility and  reversibility.  The  results  are  shown  in  Figs.  26  and 
27,  and  indicated  the  following  conclusions: 


a.  The  near-linearity  of  the  thermal  increase  of  the  elastic 
strain  justified  use  of  the  normalized  form  of  eq.  (9)  for  the 
thermoelastic  case  (Fig.  26)  . 


b.  The  measured  strain  distribution  along  the  external  EAL 
surface  was  in  fair  agreement  with  eq.  (9).  The  near-uniformity 
of  the  strain  distribution  through  the  major  midsection  was  evi- 
dent in  Fig.  27. 


The  experimental  results  provided,  however,  only  an  indirect 
and  partial  confirmation  of  the  interlaminar  shear  and  normal 
stress  distributions;  it  did  not  indicate  the  actual  interlaminar 
characteristics  close  to  the  edge,  which  is  one  of  the  main 
objectives  of  the  present  investigation. 
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5.2  Thermoelastic  Behavior  of  a Nonsymmetrical  Doubler  (NMD) 


The  following  experimental  and  analytical  work  deals  with  a 
nonsymmetrical  model  composed  of  two  strips,  each  of  a different 
material,  bonded  together  by  a thin  layer  of  polymeric  material 
and  subjected  to  thermal  changes  (see  Fig.  28) . These  changes 
and  the  different  coefficients  of  expansion,  bring  about  thermal 
stresses  causing  the  strips  to  deflect.  The  deflection  can  be 
measured  and  compared  with  analytical  solutions. 

5.2.1  Analysis. 

Equilibrium  conditions  of  the  NMD  model  give  (see  Fig.  28) : 


Nj  = N2  = N 


Mi  + M2  = N[jp-  + h0  + J*-]  = N [H] 


Vi  + V2  + xho  — 0 


This  leads  to 


The  six  equations,  together  with  the  compativility  condi- 
tions yield: 


£ i - £ 2 _ 1^  dl_ 

ho  G dx 


p being  the  beam  curvature. 


«.  i,  -ri*  * • - ■* 
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Introducing  elastic  stress-strain  relations  and  linear 
thermal  efficients  of  expansion  results  in 

(20) 

T being  the 
thermal  change 

(21) 


This  leads  to  the  following  governing  differential  equation 


[D2  - C21]M2  = C 3 1 

where 

C. 

cij  - c T- 

C‘  ' c^HG 

r = _6 6y  2 (1+u)  _ 2 (1+M) 

2 Eih2  E,h2  HE  jh i HE2h2 

C 3 = (a2-a i ) T 


(22) 

(23) 

(24) 

(25) 

(26) 


where  2c  = JL-  is  the  NMD  length  and 


y 


Ej  hj_ 
E,  h 2 


(27) 


Solving  these  governing  equations  with  the  boundary  condition 
M = M2  = 0 at  x = x/c  = ± 1 yields: 


(28) 


(29) 


Eqs.  15-17  lead  to  the  following  expressions: 


T = C3  2 (1+y) 
/CiC2  cH 


(30) 


a 


Cj 

Ci 


(u-l) 

~2c^~ 


coshVeTTx. 

cosh/C2 i 


(31) 


which  are  the  interlaminar  tensile  and  shear  stress  distribu- 
tions through  the  adhesive  layer. 

5.2.2  Discussions. 


Once  the  stress  distributions  of  t and  o are  known,  signifi- 
cant conclusions  may  be  drawn  as  to  the  effect  of  the  material 
and  geometrical  parameters  on  the  1AL  shear  and  normal  stresses 
while  subjected  to  temperature  changes. 

For  the  special  case  of  hi  = h2 , one  obtains 


c>  - (1+u> 

(32) 

„ _ 4 [12y  + 1-y 2 ] 4y 

2 E j H 2 E i H2 

(33) 

C3  = AotT 

(34) 

Assuming  = R0  to  be  the 

leads  to  the  following  equations: 

? ■ | /£±P  /± 

2 y 


rigidity  parameter  of  the  IAL 


sinh/C; i x 
cosh/C2 i 


~ _ £ - Arv  , (U-l)'  1 cosh/c77x 

u rp  V a /l  i_ , . \ / 


4 (1+y)  R cosh/c77 


Both  t and  a are  affected  by  three  main  parameters: 


(36) 
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a.  Aa,  which  is  the  measure  of  thermal  incompatibility; 

b.  a moduli  incompatibility; 

c.  the  rigidity  parameter  of  the  adhesive  layer. 

In  Figs.  29  and  30,  both  t and  a were  drawn  as  functions 
of  three  different  adhesive  rigidities  for  strips  composed  of 
aluminum  2024/T-3  and  unidirectional  (GRP)  prepreg  1002  (a  3M 
product).  (See  Table  3.) 

5.2.3  Experimental  Procedure. 

The  main  objective  of  the  experimental  phase  was  to  examine 
the  validity  of  the  analytical  solution  postulated  in  the  previous 
chapter.  For  this  purpose,  the  relation  between  the  moments  in 
the  different  strips  and  the  central  deflection  w of  a simply  sup- 
ported beam  was  calculated  and  found  to  be: 


w = 

EIlC2 


C i o2 (l-cosh/c7i 
C2  cosh/eTT 


*) 


(37) 


Special  grips  were  designed  to  enable  the  measuring  of  the  deflec- 
tion. The  testing  device,  together  witn  the  test  SDecimen,  was 
spaced  in  an  oven,  and  the  temperature  was  changed  by  10  degree 
intervals.  The  deflections  were  measured  by  a catatometer  at 
several  points  along  the  beam  (see  Fig.  31) . The  mid-point  de- 
flection was  drawn  as  a function  of  temperature.  In  Fig.  32, 
the  experimental  results  are  compared  with  the  analytical  ones; 
as  shown,  there  was  a good  agreement  between  the  two. 

5.2.4  Conclusions . 

a.  The  nonsymmetncal  multi-material  doubler  system  proved 
to  be  a good  tool  for  the  analysis  of  stress  distribution  through 
an  adhesive  layer. 

b.  The  stresses  within  the  IAL  were  reflected  in  a measure- 
able  deflection  of  the  doubler.  This  fact  may  serve  as  a simple 
yet  effective  indirect  method  for  deriving  experimentally  the 
distribution  and  critical  values  of  interlaminar  shear  and  normal 
stresses  caused  by  thermal  changes. 

c.  As  seen  from  Figs.  29  and  30,  adhesive  rigidity  parameters 
had  a major  influence  on  both  stress  peak  and  its  gradient. 

Thick  layers  of  flexible  adhesive  will  result  in  low  stresses  at 
the  edges.  On  the  other  hand,  thin  layers  of  rigid  adhesive  will 
cause  a very  high  stress  concentration  at  the  doubler  edges. 

These  stresses  are  critical  as  they  cause  premature  failure. 
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6.  SHEAR  STRESS-STRAIN  CHARACTERISTICS  OF  THE  ADHESIVE  LAYER 
IN  SITU 


6 • 1 Introduction 


The  main  problem  in  characterizing  an  adhesive  as  a structural 
material  is  whether  its  behavior  in  a confined  state  differs 
from  that  of  its  bulk,  material  reference.  To  resolve  this  prob- 
lem, which  is  essential  for  an  evaluation  of  design  allowables  of 
bonded  joints,  it  is  necessary  to  carry  out  elaborate  and  precise 
measurements  and  to  develop  a special  testing  device. 

A crucial  problem  in  such  testing  stems  from  the  fact  that 
the  thickness  of  an  adhesive  has  a major  effect  on  its  mechanical 
performance  because  of  the  adherend-adhesive  boundary  constraints 
and  the  small  measurable  displacements  obtained  in  the  case  of 
commonly  used  adhesives,  which  are  relatively  very  thin.  These 
restrictions  limited  the  available  data  [23,37,46,54]  which  origi- 
nated in  the  work  of  Hughes  and  Rutherford  [24]  . The  available 
information,  which  is  mainly  concerned  with  adhesive  behavior 
under  torsion  and  uniaxial  tension,  indicates  the  following: 
beyond  approximately  0.25  mm,  adhesive  strength  and  effective 
modulus  decrease  moderately  with  adhesive  thickness  and  tends  to 
level  off  at  a value  close  to  that  of  the  bulk  material.  Below 
this  thickness,  a steep  increase  in  adhesive  strength,  and  to  a 
lesser  degree,  of  effective  tensile  modulus,  is  evident. 

During  the  past  research  year,  shear  strength  data  of  adhesives 
in  situ  were  evaluated  by  a specially  constructed  torsion  device 
[65] . During  the  present  year,  this  device  was  improved  and  modi- 
fied by  the  attachment  of  a torsional  displacement  accessory 
(Fig.  33)  which  enabled  the  recording  of  the  shear  moment-displace- 
ment relationship  by  means  of  an  Instron  extensometer  and  recorder. 
The  testing  procedure  is  also  more  efficient  because  of  the  separa- 
tion between  the  adhesive-adherend  unit  (which  is  exchangeable)  and 
the  measuring  system. 

6.2  Testing  Procedure 

The  main  objective  of  this  test  series  was  the  evaluation 
of  the  initial  shear  characteristics  at  the  linear  elastic  range 
far  below  the  ultimate  level.  Hence,  aluminum  surfaces  were  not 
treated  but  only  cleaned  by  MEK.  The  adhesive  was  composed  of 
Shell  resin  828  and  General  Mills  Versamid  V-140,  in  a ratio  of 
2:1.  The  adhesive  was  cast  between  the  adherends'  controlled 
spacing,  which  varied  from  0.1  to  2.0  mm.  Curing  was  for  24  hours 
at  room  temperature  followed  by  48  hours  at  50 °C. 
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After  curing,  the  shear  specimens  were  installed  in  the 
torsional  device  which  was  mounted  on  the  Instron  tester.  Cross- 
head speed  was  0.5  mm/min.  For  calibration,  shear  displacement 
of  specimens  composed  of  tightened  adherends  of  zero  adhesive 
thickness,  were  recorded  to  be  deducted  from  the  over-all  shear 
displacement  of  the  bonded  specimens. 

6.3  Results  and  Discussion. 

Shear  strains  and  stresses  were  computed  from  the  recorded 
torsional  moment-displacement  curves  by  assuming  linear  shear- 
strain  distribution  throughout  the  adhesive  thickness  and  uni- 
form shear  strain  distribution  at  the  corss-sectional  area. 

Typical  shear  stress-strain  relationships  for  different  adhesive 
thicknesses  are  shown  in  Fig.  34.  Low  ultimate  values  were  ob- 
tained due  to  the  deliberately  weak  inter facial  bonding.  However, 
a qualitative  difference  can  be  distinguished  between  the  be- 
havior of  relatively  thin  and  thick  adhesives,  namely:  larger 
ultimate  strains  at  lower  stress  levels  which  may  indicate  higher 
ductility  and  are  typical  of  thin  adhesives  compared  with  their 
thick  counterparts. 

Initial  effective  shear  modulus  values,  which  were  derived 
from  the  stress-strain  curves  with  the  same  assumptions,  are 
shown  in  Fig.  35  as  a function  of  adhesive  thickness. 

In  spite  of  the  reasonable  scatter  of  the  moduli  data,  the 
general  trend  is  clear,  i.e.,  an  almost  linear  increase  of  modulus 
with  an  increase  in  thickness.  The  higher  values  of  shear  modulus, 
within  the  range  of  90-110  kg/mm2  found  for  thick  adhesive  layers 
(above  h =1.0mm)  , are  compatible  with  shear  data  obtained  for  the 
bulk  epoxy  material. 

The  conflicting  evidence  of  higher  effective  tensile  modulus 
of  thin  adhesive  layers  reported  earlier  can  be  settled  by  two 
arguments : 

a.  This  trend  was  found  to  prevail  below  0.1  mm. 

b.  The  case  of  tensile  loading  where  the  lateral  restraint 
of  metal  adherends  reduced  longitudinal  deformations,  is 
basically  different  from  the  relatively  free  and  pure 
shear  situation  which  prevailed  in  the  torsional  shear 
test. 

6.4  Future  Program. 

The  intension  is  to  continue  the  present  investigation  along 
the  following  lines: 


a.  Study  of  the  effect  of  adhesive  thickness  on  tensile 
modulus,  both  experimentally  afid  analytically,  which  will  take 
into  account  the  lateral  constraint  of  adherends  and  the  differ- 
ent characterstics  of  the  interfacial  adhesive  boundary  layers. 

b.  Taking  the  same  approach,  an  experimental  and  analytical 
study  of  the  effect  of  adhesive  layer  thickness  on  strength  in 
shear  and  uniaxial  tension 

c.  A device  for  combined  shear-axial  loading  of  adhesive 

in  situ  had  been  designed  and  will  be  utilized  for  deriving  adhe- 
sive strength  characteristics  under  combined  states  of  stress. 
Resultant  data,  it  is  hoped,  will  provide  the  basis  for  evaluat- 
ing failure  envelopes  of  adhesives  which  is  essential  for  strength 
prediction  of  the  SMD  model  ^epificaUy , and 'f or  different  adhssivs 
bonded  joint  models  in  general. 


7.  SUMMARY  AND  GENERAL  CONCLUSIONS 


The  present  report  consists  of  several  topics  concerned  with 
stress  and  strain  analysis  of  the  interlaminar  adhesive  layer 
within  a doubler  model.  The  following  conclusions  were  drawn: 

a.  Two  dimensional  stress  distribution  through  IAL  derived 
by  the  finite  element  method  revealed  the  following  character- 
istics : 

i)  Almost  uniform  shear  and  lateral  normal  stress  dis- 
tribution prevailed  throughout  the  IAL  thickness 
along  its  "middle  zone". 

ii)  Extreme  non-uniform  distributions  exist  at  the 
"boundary  zone". 

iii)  Shear  stresses  attained  their  peak  at  the  "boundary 
zone"  and  dropped  steeply  towards  zero  at  the  IAL 
edges,  which  is  consistent  with  boundary  and 
equilibrium  conditions. 

iv)  Axial  normal  stresses,  which  were  negligible  at  the 
IAL  "boundary  zone",  attained  the  same  order  of 
magnitude  as  other  interlaminar  stresses  along  the 
"middle  zone"  close  to  the  central  adherend. 

v)  The  principal  interlaminar  stresses  reached  their 
extreme  value  close  to  the  intersection  of  the  IAL 
edge  and  the  IAL  central  adherend  interface. 


b.  Stress  distributions  along  the  IAL  which  were  determined 
by  a direct  measurement  of  shear  and  normal  relative  displacements, 
show  good  agreement  with  both  analytical  closed  form  solution  and 
with  numerical  solution  derived  by  finite  element  method. 

c.  Thermoelastic  investigation  of  IAL  stresses  within  sym- 
metrical and  nonsymmetrical  doubler  models  showed  good  agreement 
between  experimental  measurements  and  analytical  closed  form 
solutions . 

d.  The  stress-strain  relationship  of  the  adhesive  layer 
in  situ,  when  measured  by  a special  torsion  device,  indicated 
that  the  adhesive  shear  modulus  increased  moderately  with  the 
thickness  of  the  adhesive  and  approached  the  respective  modulus 
of  the  bulk  epoxy  at  a thickness  of  about  1.0  mm. 


8.  RECOMMENDATIONS  FOR  FUTURE  RESEARCH 


It  is  recommended  that  the  present  research  be  continued 
along  the  following  lines  (see  Figs.  36,37): 

a.  Parametric  study  of  the  effects  of  variations  in  IAL 
thickness  and  its  edge  geometry,  as  well  as  its  material 
characteristics  (by  FEM) . 

b.  Two  dimensional  stress  analysis  of  SMD  at  the  nonlinear 
range  of  IAL  behavior  (by  FEM) . 

c.  Study  of  the  elastic  and  thermoelastic  three  dimensional 
behavior  of  SMD. 

d.  Investigation  of  strength  and  failure  characteristics  of 
IAL  within  SMD  as  related  to  respective  behavior  of 
adhesive  in  sit u,  under  combined  loading. 

e.  Investigation  of  the  environmental  effects  (temperature 
and  humidity)  on  IAL  mechanical  behavior,  as  reflected 
by  SMD  performance. 
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GLOSSARY 

C. . - Structural  parameters  of  NMD 

c - lenglth  of  IAL  (c  = i/2) 

Eo ,Ei ,E2 , (E^)  - Young's  moduli  of  IAL,  CAL  and  EAL,  respectively 

F - axial  applied  force  per  unit  width 

Go  - shear  modulus  of  IAL 

H - measure  of  NMD  thickness 

ho ,hi ,h2  - thickness  of  IAL,  CAL  and  EAL,  respectively 

K.  - constants 

1 

i - IAL  length  (&= 2c) 

M2  =M  - moment  function  along  the  EAL 

N i , N 2 =N  - normal  force  along  the  CAL  and  EAL,  respectively 

Pc  - effective  applied  axial  normal  stress 

po  - reference  axial  stress  acting  on  EAL  under  uniform 

axial  strain 

P - axial  applied  force 

R - structural  parameter  of  NMD 

T - temperature  difference 

u - displacement  in  x-direction 

V 2 =V  - shear  force  along  EAL 

w - displacement  in  z-direction 

x,y,z  - axial,  transverse  and  lateral  coordinates,  respectively 
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- coefficient  of  thermal  expansion 
Y0  - shear  strain  at  1AL 

6 - material  and  geometrical  constant 

eo,Ei,e2  - axial  normal  strain  at  IAL , CAL  and  EAL,  respectively 

ri  - material  and  geometrical  constant 

X - material  and  geometrical  constant 

Vi  - material  and  geometrical  constant 

v,vo  - Poisson's  ratio  of  adherend  and  adhesive  layers  respectively 

p - curvature  of  beam 

Tzxo=To  " ax;'-a-'-  shear  stress  (in  x-direction)  activing  at  IAL 

T - transverse  shear  stress  (in  y-direction)  acting  at  IAL 

zy° 

oz0  =Oo  - lateral  normal  stress  (in  z-direction)  acting  at  IAL 

a , a ,cr  - longitudinal  normal  stresses  acting  at  IAL,  CAL 

1 2 and  EAL  respectively 


yo 

X 

SMD 

NMD 

IAL 

CAL 

EAL 

FRP 

SLJ 

DLJ 

FEM 


a.. 


yi  yz 


- transverse  normal  stresses  acting  at  IAL,  CAL 


and  EAL  respectively 

- material  and  geometrical  constant 

- non-dimensional  axial  coordinate 

- symmetrical  doubler  model 

- nonsymmetrical  doubler  model 

- interlaminar  adhesive  layer 

- central  adherend  layer 

- external  adherend  layer 

- fiber  reinforced  plastics 

- single  lap  joint 

- double  lap  joint 

- finite  element  method 
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FIG.  8 AXIAL  DISTRIBUTION  OF  SHEAR  STRESSES  THROUGH 
THE  EXTERNAL  ADHERENDS . 


finite  element  1 


NON  DIMENSIONAL  AXIAL  COORDINATE 
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FIG.  12  AXIAL  DISTRIBUTION  OF  SHEAR  STRESSES  THROUGH  THE 


0.20 


47 


NORMALIZED  INTERLAMINAR  AXIAL  NORMAL  STRESS 


NONCXMENSIONAL  AXIAL  COORDINATE 


NORMALIZED  INTERLAMINAR  SHEAR  STRESS  T«0 


NORMALIZED  INTERLAMINAR  PRINCIPAL  STRESS  3j, 
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FIG.  20  ELECTROMECHANICAL  EXTENSOMETER  FIG.  21  ELECTROMECHANICAL  EXTENSOMETER 

FOR  MEASURING  SHEAR  DISPLACEMENT.  FOR  MEASURING  LATERAL  NORMAL 

DISPLACEMENT. 


experimental  (specimen  no.  A21/1) 


FIG.  23  INTERLAMINAR  SHEAR  STRESS  DISTRIBUTION  (T 
EXPERIMENTAL  VS  ANALYTICAL  SOLUTION 
(MODEL  A21) . 


experimental  (specimen  no.  A21/1) 


INTERLAMINAR  LATERAL  NORMAL  STRESS  DISTRIBUTION  (a 
EXPERIMENTAL  VS  ANALYTICAL  SOLUTION  (MODEL  A21) . Z° 


FIG.  27  NORMALIZED  STRAIN  DISTRIBUTION  RESULTING  FROM  THERMAL 
CHANGES  ALONG  THE  EAL  EXTERNAL  SURFACE  (MODEL  A02) . 


FIG.  28  STRESS  ANALYSIS  SCHEME  OF  A NONSYMMETRICAL 
DOUBLER  MODEL  (NMD) . 


FIG.  31  SCHEMA  FOR  MEASURING  SYSTEM  OF  NMD  DEFLECTIONS 


FIG.  32  THERMOELASTIC  DEFLECTION  CURVE  OF  NMD 


FIG.  33  SCHEMA  OF  MODIFIED  TORSION  TESTING  DEVICE  FOR  DIRECT 
MEASUREMENT  OF  THE  SHEAR  STRESS -STRAIN  BEHAVIOR  OF 
THE  ADHESIVE  LAYER  IN  SITU. 


TORSIONAL  SHEAR  STRESS  T (kgAnm*) 


THICKNESS  OF  ADHESIVE  LATER  h0  (mm) 

FIG.  35  THE  EFFECT  OF  ADHESIVE  LAYER  THICKNESS  ON 
INITIAL  SHEAR  MODULUS. 


FIG.  36  FLOW  CHART  FOR  A FUTURE  PROGRAM  ON  NON-LINEAR  BEHAVIOR  OF  THE  INTERLAMINAR 
ADHESIVE  LAYER  WITHIN  A MULTIMATERIAL  DOUBLER  MODEL. 


FIG.  37  FLOW  CHART  FOR  A FUTURE  PROGRAM  ON  THREE  DIMENSIONAL  BEHAVIOR  OF  THE  INTERLAMINAR 
ADHESIVE  LAYER  WITHIN  A MULTIMATERIAL  DOUBLER  MODEL. 


TABLE  1:  DIMENSIONS  AND  ELASTIC  PARAMETERS  OF  TEST  MODELS. 

Dimensions  [mm] 
io  hi  h2  c 

1.0  1.5  1.0  65 

1.0  1.5  1.0  65 


Elastic 

Moduli 

[kg/mm2 ] 

Go 

E0 

Ei 

e2 

90 

i 250 

7500 

7500 

90 

' 250 

3500 

7500 

TABLE  2:  ELASTIC  AND  THERMOELASTIC  PARAMETERS  OP  TEST  MODEL  AO  2 


Direction 


COMPOSITION  OF  SMD  (TABLES  1,2) 


Adhesive 


2:1  mix  of  epon  828 
and  versamid  140 


TABLE  3:  DIMENSIONS  AND  ELASTIC  PARAMETERS  OF  NMD  TEST  MODEL. 


Elastic 

Moduli 

kg/mm2 

Dimensions  (mm) 

f ~ 

G0 

Eo 

I t 

Ei 

e2 

ho 

hi  h2  c 

90 

250 

7500 

3500 

0.1 

1 1 90 

i 

Composition; 
Adherend  1: 
Adherend  2; 

Adhesive: 


Aluminum  2024/T3  non-clad. 

E glass  fiber  reinforced  unidirectional 
laminate . 

2:1  mix  of  Epon  828  resin  and  Versamid 
140  hardener. 
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FINITE  ELEMENT  COMPUTER  PROGRAM  FOR  ELASTIC  STRESS 

analysis  in  orthotropic  bonded  structural  systems. 
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